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Abstract 



The effect of surface fluctuation on the observables of quark-hadron phase 
transition is studied. The Ginzburg-Landau formalism is extended by the in- 
clusion of an extra term in the free energy that depends on the vertical dis- 
placements from a flat surface. The probability that a bin has a particular net 
. displacement is determined by lattice simulation, where the physics input is 

color confinement. The surface fluctuation from bin to bin is related to multi- 
plicity fluctuation, which in turn is measured by the factorial moments. It is 
. found that both the F-scaling behavior and the scaling exponent are essentially 

unaffected by the inclusion of surface fluctuation. 



1 Introduction 
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Scaling behaviors of hadronic observables in quark-hadron transition (PT) have been 
studied in the framework of Ginzburg-Landau (GL) theory [0-0]. Since it is a mean- 
field theory, it cannot account for the dynamical fluctuations that can occur in heavy- 
ion collisions. An attempt to amend that defect has been made by treating the prob- 
lem in the Ising model, so that by lattice simulation the effect of spatial fluctuation 
of hadron density at the surface of quark-gluon plasma can be considered || . In this 
paper we treat the problem of spatial fluctuation in a very different way. We gener- 
alize the GL description of second-order PT in 2D by including a term that accounts 
for surface fluctuation. The deviation from a flat surface is simulated by use of a 
term in the free energy that represents the confining potential between quarks. By 
associating bulges on the surface to regions where hadronization can take place more 
readily, we can determine the effect of spatial fluctuations on the hadronic observ- 
ables. Thus this work represents a significant step toward making the GL formalism 
a more realistic description of the PT problem in heavy-ion collisions. 

In Refs. [DJ-pj it has been found that the scaling behavior characteristic of the 
PT is associated with the normalized factorial moments F q . Actually, that scaling 
behavior does not exhibit intermittency |J, which is the power law F q oc M Vq , where 
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M is the number of bins in a fixed volume of phase space. Instead of M scaling, 
which we use to refer to the M"- Pq behavior, we have found that pure GL leads to F 
scaling, which is 

F q oc F%\ (1) 
when M is varied. Furthermore, j3 q satisfies the simple formula [|TJ 

P q ={q-l) v . (2) 

It is this v that has been referred to as the scaling exponent. It is not to be confused 
with the conventional critical indices in statistical physics, to which v bears absolutely 
no relationship. In problems where the temperature T is not observable, but the 
multiplicities of particles produced are (such as photons in lasers |7] and pions in 
heavy- ion collisions |4[]), F scaling has been found to be the only scaling behavior 
that can be put to experimental tests. Thus far the value of u at 1.3 has proven 
to be remarkably universal theoretically, and has been verified experimentally to be 
accurate in quantum optics @. In this paper we shall check whether F scaling remains 
valid when spatial fluctuations are taken into account, and if so, whether v = 1.3 is 
still true. 



2 Formulation of the problem 

We give first a very brief review of the usual method in the GL approach to hadronic 
observables The multiplicity distribution P n is expressed as a functional inte- 

gral 

Pn = y V<f>± (/ rf 2 *|0| 2 ) V/^'V^ (3) 

where Z = J T>cf) e~ F , and the GL free energy is 

F = J d 2 x(a\<P\ 2 + b\<P\ A + c\V<P\ 2 ). (4) 

Here the spatial integration is over a 2D space, which may be taken to be the surface 
of the cylinder containing quark-gluon plasma. When the surface temperature T, on 
which the parameters a, b, and c in (f|) depend, is low enough, hadrons are formed on 
that surface and are removed from the confining medium that exists within. If the 
integration in x is limited to a small bin of size A, then the A-dependence of P n can 
best be determined in terms of F q , where 

_ (7i(w-l)---(w-g+l)) 
9 ~ {n) q 1 ' 

the average being defined by use of the distribution P n . It is this F q that has the 
scaling properties ([I]) and (0). 
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In most treatments PJ-H the third term in (H) is neglected. Its inclusion is 
considered by Hwa and Pan [||, and is found to lead to only a small change in the 
value of v. Our approach here is to replace that term by another that has explicit 
spatial dependence with dynamical origin. But before embarking on that, let us 
rewrite the formalism in a more compact way. 

If we write (|5|) in the form 

F q = fjfl (6) 

then we have 

h = E 7 '^ p n = / dm m«D q (m), (7) 

^ q {n-q)\ J 

where D q {m) is proportional to e~ F ^ m \ and 

m = [ d 2 x\<f)(x)\ 2 = A\<f)\ 2 . (8) 

Without the last term of (|), may be regarded as constant inside A so the last 
expression in (§) follows. With the interpretation of |0| 2 as the hadron density, an 
order parameter, we may regard m as the "dynamical" multiplicity in a bin of area 
A. The statistical fluctuation represented by the Poissonian factors in front of e~ F 
in (Q) is removed by taking the factorial moment, so that the series in terms of the 
"experimental" multiplicity n in ([?]) is replaced by the integral over the "dynamical" 
multiplicity m, which need not be an integer. f q has multiplicative factors that are 
cancelled by the same factors in ff , when the ratio in (|6]) is taken. Leaving out those 
factors, we can write f q more transparently as 

POO 

f q = / dmm q exp[-(m-m a ) 2 /(A/b)], (9) 
Jo 

where 

m a = -Aa/2b. (10) 

PT occurs when a < 0, so m a is always positive in the following consideration. Note 
that a is dimensionless, while b has the dimension of area. Thus the physical A always 
appears in the ratio A/b, i.e., in the GL description the area is measured in units of 
b. In [0] it is shown that F q is a function of only one composite variable, Aa 2 /2b. 
In the following because of our way of introducing spatial fluctuations, there will be 
seperate dependences on a and A/b. Hence, @ will be our starting point, when there 
is no spatial fluctuation. 

To introduce surface fluctuations, let us first map the 2D surface to a square 
lattice of size L 2 , on which we impose periodic boundary condition. With bins of size 
5 2 , there are M = (L/5) 2 bins on the lattice. At each lattice site we assign a variable 
Zi that represents the vertical displacement from the flat surface at site i. We restrict 
Zi to be in the range — 1 < Zi < 1, deferring the physical scale of perpendicular 
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displacement to another parameter to be introduced below. Zero displacement is 
defined by the global constraint 

where the sum is over all sites on the lattice. 

The physics governing the ^-displacement is color confinement. In the discretiza- 
tion of the physical space we do not mean that the lattice spacing is to be identified 
with the actual average distance between partons on the surface of the plasma. The 
lattice is a coarse-grain representation of the surface, and we let each site carry the 
net color charge of the basic cell (of area lattice- spacing squared) surrounding the site. 
For brevity we call such a color charge at the site a quark. Our aim is to simulate 
vertical fluctuations of the surface. With the lattice spacing being fixed, the vertical 
displacement of a site relative to its neighbors represents further separation among 
the quarks. Since energy is required to increase the quark separation, we add a term 
to the free energy to inhibit vertical displacement, calling it F s to represent spatial 
contribution, 

F s = Cj2\zi-Zj\, (12) 

m 

where the sum is over the nearest neighbors j of each site i, and then over all i. 
The confinement parameter C includes the (fcgT) -1 factor, as in (ED), and the scale 
factor for perpendicular displacement mentioned earlier. The \zi — Zj\ dependence 
represents the linearly increasing portion of the confinement potential. The — 1/r 
Coulombic portion is omitted partly because the lattice spacing imposes a minimum 
distance between quarks that prevents r from becoming small, and partly because we 
want to avoid introducing another parameter. In place of the last term of (Ejj), F s as 
given in fll^) provides a sensible description of spatial dependence of the energy of 
the system that incorporates the color dynamics into the GL formalism. A key link 
yet to be specified is the way in which Zi is related to the order parameter <fi or the 
dynamical multiplicity m. 

Our procedure is to use the Boltzmann factor exp(— F s ) to simulate a configuration 
of Zi throughout the lattice. Clearly, because of thermal fluctuation the vertical 
displacement need not be Zi = for all i, even though that is the most favored 
configuration. Because of (JTTJ) there are as many net positive displacements as there 
are negatives on the lattice. However, as we focus on a particular bin of area 6 2 , there 
may be a nonvanishing displacement for the bin 

c = E *■ (13) 

i<=6 2 

We identify a positive ( as having a net outward protrusion of the surface at the bin, 
and a negative £ as having a net inward indentation. Nucleation dynamics suggests 
that hadronization is more likely to occur at the location of protruding surface and 
less likely at where the surface is indented. Thus given a value of (f> for a bin, the 
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dynamical multiplicity should deviate from m in a way that is proportional to £. 
That is, we should consider a convolution G ® S, where G is the GL term given by 
the exponential in (0), and S represents the spatial fluctuation specified by S(Q, 
yet to be determined, which describes the probability that a bin has a net vertical 
displacement (. 

More precisely, we upgrade (|9p to include spatial fluctuations by writing 

roc 

f q = / dfifiWifi), (14) 

JO 

where the dynamical distribution D(fi) is 

/•CO 

D{n) = G®S= dm G(m)S(fj, - m) (15) 
Jo 

and 

G(m) = exp[-(m - m a f/{A/b)\. (16) 

In the next section we shall determine the distribution in (. What remains to be spec- 
ified is the relationship between ( and hadron multiplicity. Without getting involved 
with the details of nucleation, surface tension, etc., let us introduce a new parameter 
a that summarizes all the dynamical details relating vertical surface displacement to 
multiplicity fluctuation, /i — m, so that 

fi = m + a(. (17) 

Recall that m is the mean multiplicity in the GL theory. The spatial fluctuation 
of the multiplicity distribution, S(a(), to be determined in the next section is what 
appears in ( |i~5"D for the net dynamical distribution -D(^). 



3 Simulation on the lattice 

Our first task is to reproduce to GL result |l| (without surface fluctuation) by lattice 
simulation. In the analytical method used in M, there is only one composite variable, 
x = Aa 2 /2b, that F q depends on. F scaling is found to hold in the range 0.1 < lnF 2 < 
0.44, which corresponds roughly to 3 < Initio < 9. From Fig. 1 of one can see 
that such ranges of F q are for x in the range —2 < — lnx < 4. Since x is a composite 
variable, it does not matter how a, b and A are separately varied. Phenomenologically, 
a and b are not under experimental control, so A must be varied in the appropriate 
range such that the data on F q are in the range that exhibits F scaling. Now, on the 
lattice we must work with specific bin sizes 5. That necessity destroys the advantage 
of dealing with just one x variable in the analytic method. We make the identification 

A/b = 5 2 (18) 

and choose a range of S that can cover the scaling range in x mentioned above, i.e., 
0.018 < x < 7.5. The ratio of the extrema of that x range corresponds, for a held 
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constant, to the ratio of the extrema of 5 being very nearly 20 : 1. Thus, if we let the 
smallest bin to consist of just two sites in each direction, then the range of 6 should 
be 2 < 5 < 40. Choosing the lattice size to be L = 120, we have for the number of 
bins, M, on the lattice to vary from 60 2 down to 3 2 , a range that is wide enough to 
test scaling. The corresponding value of a is roughly —0.1. 

Applying the above parameters to (|9]) and (|10l), and using Monte Carlo method to 
calculate f q and then F q , we find the result in Fig. 1 shown in log-log plot. The linear 
dependences for 3 < q < 6 confirm the F-scaling behavior. Moreover, the scaling 
exponent v defined in (0) is found to be v = 1.308, as shown in Fig. 2, in essential 
agreement with the result obtained in I], as it should. 

Proceeding now to the main task of including spatial fluctuations, we start by 
setting C — 1 in ( |T2"D initially. Other values of C will be considered later. We 
use the Metropolis algorithm to simulate the configurations of z^ on each site of 
the lattice, with the Boltzmann factor e~ Fs as the weight. We discard the first 800 
sweeps for initialization. From 10 3 such configurations we calculate the distribution 
in (. However, owing to (|T7|) it is more useful to plot the distribution as a function 
of a(. The arbitrarily chosen value of a affects only the scale of the abscissa and 
the normalization of S(a(). Fig. 3 shows the results for 5 = 2 and 40, and for 
a = 1/4, which is chosen so that for the smallest bin the maximum multiplicity due 
to surface fluctuation is a( = 1. As Fig. 3(a) shows, the distribution is very nearly 
Gaussian with half-width around 0.3, plus a sharp peak at a( = 0, corresponding 
to the preponderant probability that Zi = for all four sites in a bin. On the other 
hand, for the largest bin S = 40, the distribution shown in Fig. 3(b) is negligible for 
a( > 50, although the maximum possible value, for all Zi = +1, is a( = 40 2 /4 = 400. 
Apart from the peak at a( = the half-width is around 20. The statistics is poorer in 
this case because the lattice has far fewer bins per configuration when 5 is increased 
twenty-fold. As the area 5 2 is increased by a factor of 400, the mean multiplicity 
fluctuation is increased by a factor of only 20/0.3 ~ 60. Thus the deviation from a 
flat surface does not have long wavelength modes, which would have resulted in large 
C configurations in large bins. On the other hand, the surface fluctuations are not all 
local, lest Fig. 3(b) would not be as broad as it is. 

Putting these S(a() distributions for various values of S between 2 and 40 in ( |I~5D 
and performing the convolution, we obtain the dynamical multiplicity distribution 
D(fi). Fig. 4 shows the results for a range of values of 5 with C = 1 and a = |. From 
the D(/j,) distribution we can use (|]) and (|T4|) to calculate F q , which we present in 
the logFq vs logF 2 form in Fig. 5. We see that the straightline behavior supports F 
scaling. As a is varied from 1/4 to 2, the ranges of F q values are reduced but the slopes 
of the straightlines are essentially unchanged. Thus the first part of our objective has 
been accomplished; that is, F scaling persists even when surface fluctuations are 
considered. 

Before we investigate the second part, viz., the effect on scaling exponent, we 
examine the dependence on C . The foregoing analysis is done for C = 1, a value 
that has been arbitrarily chosen. We now vary C to see the effect on S((). From 
the simulated result on S(Q for every set of values of C and a, we calculate the 
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width = (( 2 ) and show in Fig. 6 the dependence of on C for the two extreme 
bin-sizes, 5 = 2 and 40. We see that for 5 = 40 there is a maximum at C = 2, 
so the previous result on C = 1 is not far from maximum surface fluctuation. 
is lower at both higher and lower values of C because high C stiffens the surface 
and suppresses (, while low C leads to random Zi and therefore also reduces £. The 
same is not true for 5 = 2, since even in the limit of C — > the smallness of the bin 
size makes the cancellation of random Zi ineffective, resulting in a Possionian S(() 
that has nonvanishing width. As C increases, the surface becomes smoother, and Tq 
continues to decrease monotonically. To summarize, we need only consider the range 
< C < 2 for maximum effect due to surface fluctuation. 

We now can repeat the previous calculations for other values of C and produce 
figures similar to those shown in Figs. 4 and 5. However, since the results are very 
similar, we omit those figures. 

To determine the slopes (3 q of the F-scaling behavior, we fit all the points in Fig. 
5 for a given q value, including all values of a, by a straight line. The resultant j3 q 
is then shown in Fig. 7 for a particular value of C {C = 1 for Fig. 5). The error 
bars in Fig. 7 correspond to the deviations from the straightline fits of the points in 
Fig. 5 for C = 1 only. For other values of C the results from similar calculations are 
also shown in Fig. 6. No corresponding error bars are shown, since they are similar. 
The f3 q values are determined only for integer values of q, but since the dependence 
on C is so insignificant, we have for clarity exhibited the j3 q values by straightline fits 
according to (fj). Even so, the overlap obsures the individual lines for different values 
of C, but successfully demonstrates that the scaling exponent v is insensitive to C. 
The net value of v that can be concluded from Fig. 7 is 

v = 1.306 ±0.035, (19) 

when surface fluctuations are taken into account. Compared to the result of v = 1.308 
obtained from Fig. 2 for pure GL, we see that the change in v is well within the errors. 

4 Conclusion 

To study the effect of surface fluctuation on the observables of quark-hadron PT, 
we have incorporated into the GL description an additional term in the free energy 
that depends on the vertical displacements from a flat surface. We have used lattice 
simulation to determine the probability that a bin has a net displacement (. The 
physics that controls such displacements is color confinement. We connect the surface 
fluctuations to multiplicity fluctuations, which in turn are quantified by normalized 
factorial moments F q . The F scaling of F q , when the bin size is varied, results in the 
determination of the scaling exponent v. The aim of this work has been to study the 
dependence of v on surface fluctuations, after verifying that the scaling behavior is 
not destroyed. 

There are two parameters in the calculation, on which we have no a priori informa- 
tion. One is C; the other is a. C depends on the strength of confinement, and relates 
the vertical displacement Zi (arbitrarily normalized) to the free energy (normalized by 
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the transition temperature). The parameter a relates the net displacement ( of a cell 
to the deviation of the cell multiplicity from the mean. We have studied the variation 
of u, when C and a are varied over the important ranges < C < 2 and < a < 2. 
When C ~ 2 or a ~ 2, we do see significant fluctuations. Our result, however, shows 
that z/ is essentially unchanged when C and a are varied in the ranges considered. 
There is no doubt that when a is large enough, the scaling behavior will break down. 
But then our method of calculation should become invalid in that case, since large 
multiplicity fluctuations cannot be treated in the framework of the mean-field theory. 
Besides, there is no physical reason to expect that small surface fluctuation can lead 
to large multiplicity fluctuations. 

The significance of this work is therefore in finding that the scaling behavior of 
pure GL is stable against the perturbation introduced by surface fluctuations. Thus 
this result improves the Ginzburg-Landau description of quark-hadron phase transi- 
tion and provides a more physical connection to the realistic situation in heavy-ion 
collisions. The end result is that F scaling persists even when surface fluctuations 
are taken into account, and that the scaling exponent, which is experimentally mea- 
surable, remains at v ~ 1.3. 
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Figure Captions 

Fig. 1 F scaling in Ginzberg-Landau theory without spatial inhomogeneity, i.e., 
c = in Eq. (4). 

Fig. 2 Log- log plot verifying Eq. (2), using the data of Fig. I. 

Fig. 3 Distribution of surface fluctuation for C — 1 and a = 1/4. The bin sizes 
are different in (a) and (b) as indicated. 

Fig. 4 Distribution of dynamical multiplicity when the GL description of PT is 
supplemented by surface fluctuations. 

Fig. 5 F-scaling plot of F q for C = 1 and a range of a. 

Fig. 6 Width of S(() vs C for 5 = 2 (right scale) and 40 (left scale). 

Fig. 7 The slopes f3 q for a range of values of C. Although (3 q are determined for 
integer values of q, they are plotted as straight lines to avoid overlap of points. 
The error bars are for C — 1 only, the ones for other values of C being similar. 
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Fig 5. 
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